C SOURCES:
C    Thomas  M.  Henderson,  Benjamin  G.  Janesko,  and  Gustavo  E.  Scuseria.
C    Generalized  gradient approximation model exchange holes for range-separated hybrids.
C    The Journal of Chemical Physics,128(19):194105, may 2008.



C*****************************************************************************
      pure subroutine ESRX_wPBE(rho_a, gamma_aa, mu, Ea)
C*****************************************************************************
C   Implemented by E.R. Kjellgren.
C
C   Subroutine generated using Sympy 1.3
C   Generated: March 22, 2019
C*****************************************************************************
      implicit none
      real*8, intent(in) :: rho_a, gamma_aa, mu
      real*8, intent(out) :: Ea
      real*8 :: x0, x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x
     &13, x14, x15, x16, x17, x18, x19, x20, x21, x22, x23, x24, x25, x2
     &6, x27, x28, x29, x30, x31, x32, x33, x34, x35, x36, x37
      Ea = 0.0d0
      x0 = rho_a**1.33333333333333d0
      x1 = 0.217336917462899d0
      x2 = gamma_aa*x1
      x3 = rho_a**(-2.66666666666667d0)
      x4 = x2*x3
      x5 = 0.00478594012208448d0*gamma_aa**3.5d0*rho_a**(-9.333333333333
     &33d0)
      x6 = 0.0102659822546843d0*gamma_aa**3*rho_a**(-8.0d0)
      x7 = 0.0220208337265180d0
      x8 = gamma_aa**2.5d0*rho_a**(-6.66666666666667d0)*x7
      x9 = 0.0472353356922751d0*gamma_aa**2*rho_a**(-5.33333333333333d0)
     &
      x10 = 0.101321183642338d0
      x11 = gamma_aa**1.5d0*rho_a**(-4.0d0)*x10
      x12 = 0.466194077035412d0
      x13 = sqrt(gamma_aa)*x12
      x14 = x3*(0.00177707291666667d0*x11 + 0.00121096701974236d0*x4 + 5
     &.04056248768705d-6*x5 - 4.21555121527778d-5*x6 + 0.000240776365358
     &387d0*x8 - 0.000919311578213079d0*x9)/( -5.64967128377635d-9*gamma
     &_aa**4*rho_a**(-10.6666666666667d0) + 3.14930255098365d-10*gamma_a
     &a**4.5d0*rho_a**(-12.0d0) + 0.231225d0*x11 - 0.944751281556643d0*x
     &4 + 3.67304772509327d-5*x5 - 0.000264921875d0*x6 + 0.0027046704323
     &9359d0*x8 - 0.0292938845354061d0*x9 + 1.0d0 + 1.46748378204436d0*x
     &13/x0)
      x15 = 0.0757133580346725d0*x14
      x16 = x15*x2
      x17 = x16 + 0.60965d0
      x18 = 0.30285343213869d0*mu**2*rho_a**(-0.666666666666667d0)
      x19 = x1*x18
      x20 = x17 + x19
      x21 = sqrt(x20)
      x22 = mu*rho_a**(-0.333333333333333d0)
      x23 = x12*x22
      x24 = x23/x21
      x25 = sqrt(gamma_aa*x15 + x18)
      x26 = x16 + 0.757211d0
      x27 = sqrt(x19 + x26)
      x28 = 0.550321208149104d0*x22
      x29 = x12*x28
      x30 = 1d0/(x21 + x29)
      x31 = x14*x2
      x32 = 0.151426716069345d0*x31
      x33 = x17**2
      x34 = x4/(0.0189283395086681d0*x4 + 1.0d0)
      x35 = mu**3*1d0/rho_a*x10*x20**(-1.5d0)
      x36 = x17**3
      x37 = x17**3.5d0
      Ea = -1.8610514726982d0*x0*(1.10064241629821d0*x23*(x12*x25 - x27)
     & + x32*log( x12*x30*(x25 + x28)) - (x32 + 1.514422d0)*log(x30*(x27
     & + x29)) + 0.757211d0 - (0.0189283395086681d0*mu**5*rho_a**( -1.66
     &666666666667d0)*x20**(-2.5d0)*x7 + 1.03185226527957d0*x24 - 0.2083
     &33333333333d0*x35 - 1.0d0)*(-0.888888888888889d0*x17*( 0.015142671
     &6069345d0*x31 + 0.00112167937829144d0*x34 + 0.0474594676428533d0) 
     &- 0.0252121521774518d0*x33 + 0.807691733333333d0*x36 + 2.133333333
     &33333d0*x37*( 0.275160604074552d0*x13*sqrt(x14) - sqrt(x26)) + 1.2
     &6041162731059d0*x37)/x36 + (-0.825481812223657d0*x24 + 0.083333333
     &3333333d0*x35 + 1.0d0)*(0.0168251906743717d0*x31 + 0.0012463104203
     &2383d0*x34 + 0.0527327418253926d0)/x33 - ( 0.0260152163368769d0*x2
     &4 - 0.0472727853327222d0)/x17)
      end subroutine


C*****************************************************************************
      pure subroutine D1ESRX_wPBE(rho_a, gamma_aa, mu, Ea, d1Ea)
C*****************************************************************************
C   Implemented by E.R. Kjellgren.
C
C   Subroutine generated using Sympy 1.3
C   Generated: March 22, 2019
C*****************************************************************************
      implicit none
      real*8, intent(in) :: rho_a, gamma_aa, mu
      real*8, intent(out) :: Ea, d1Ea(4)
      real*8 :: x0, x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x
     &13, x14, x15, x16, x17, x18, x19, x20, x21, x22, x23, x24, x25, x2
     &6, x27, x28, x29, x30, x31, x32, x33, x34, x35, x36, x37, x38, x39
     &, x40, x41, x42, x43, x44, x45, x46, x47, x48, x49, x50, x51, x52,
     & x53, x54, x55, x56, x57, x58, x59, x60, x61, x62, x63, x64, x65, 
     &x66, x67, x68, x69, x70, x71, x72, x73, x74, x75, x76, x77, x78, x
     &79, x80, x81, x82, x83, x84, x85, x86, x87, x88, x89, x90, x91, x9
     &2, x93, x94, x95, x96, x97, x98, x99, x100, x101, x102, x103, x104
     &, x105, x106, x107, x108, x109, x110, x111, x112, x113, x114, x115
     &, x116, x117, x118, x119, x120, x121, x122, x123, x124, x125, x126
     &, x127, x128, x129, x130, x131, x132, x133, x134, x135, x136, x137
     &, x138, x139, x140, x141, x142, x143, x144, x145, x146, x147, x148
     &, x149, x150, x151, x152, x153, x154, x155, x156, x157, x158, x159
     &, x160, x161, x162, x163, x164, x165, x166, x167, x168, x169, x170
     &, x171, x172, x173, x174, x175, x176, x177, x178, x179, x180, x181
     &, x182, x183, x184, x185, x186, x187, x188, x189, x190, x191, x192
     &, x193, x194, x195, x196, x197, x198, x199, x200, x201, x202, x203
     &, x204, x205, x206, x207, x208, x209, x210, x211, x212, x213, x214
     &, x215, x216, x217, x218, x219, x220, x221, x222, x223, x224, x225
     &, x226, x227, x228, x229, x230, x231
      Ea = 0.0d0
      d1Ea(:) = 0.0d0
      x0 = 0.217336917462899d0
      x1 = rho_a**2.66666666666667d0
      x2 = 1d0/x1
      x3 = x0*x2
      x4 = 0.944751281556643d0*x3
      x5 = rho_a**(-12.0d0)
      x6 = 0.00104016147329585d0
      x7 = gamma_aa**4.5d0*x6
      x8 = rho_a**(-10.6666666666667d0)
      x9 = 0.00223117693796192d0
      x10 = gamma_aa**4*x9
      x11 = rho_a**(-9.33333333333333d0)
      x12 = 0.00478594012208448d0
      x13 = gamma_aa**3.5d0
      x14 = x12*x13
      x15 = x11*x14
      x16 = rho_a**(-8.0d0)
      x17 = 0.0102659822546843d0
      x18 = gamma_aa**3
      x19 = x17*x18
      x20 = x16*x19
      x21 = rho_a**(-6.66666666666667d0)
      x22 = 0.0220208337265180d0
      x23 = gamma_aa**2.5d0
      x24 = x22*x23
      x25 = x21*x24
      x26 = rho_a**(-5.33333333333333d0)
      x27 = 0.0472353356922751d0
      x28 = gamma_aa**2
      x29 = x27*x28
      x30 = x26*x29
      x31 = rho_a**(-4.0d0)
      x32 = 0.101321183642338d0
      x33 = gamma_aa**1.5d0
      x34 = x32*x33
      x35 = x31*x34
      x36 = sqrt(gamma_aa)
      x37 = 0.466194077035412d0
      x38 = rho_a**1.33333333333333d0
      x39 = 1d0/x38
      x40 = x37*x39
      x41 = -gamma_aa*x4 - 2.53214847628224d-6*x10*x8 + 3.67304772509327
     &d-5*x15 - 0.000264921875d0*x20 + 0.00270467043239359d0*x25 - 0.029
     &2938845354061d0*x30 + 0.231225d0*x35 + 1.46748378204436d0* x36*x40
     & + 3.02770543981481d-7*x5*x7 + 1.0d0
      x42 = 1d0/x41
      x43 = 0.00121096701974236d0*x3
      x44 = gamma_aa*x43 + 5.04056248768705d-6*x15 - 4.21555121527778d-5
     &*x20 + 0.000240776365358387d0*x25 - 0.000919311578213079d0*x30 + 0
     &.00177707291666667d0*x35
      x45 = x42*x44
      x46 = x2*x45
      x47 = gamma_aa*x46
      x48 = 0.0757133580346725d0*x47
      x49 = x0*x48
      x50 = x49 + 0.60965d0
      x51 = 1d0/x50
      x52 = mu**2
      x53 = 0.30285343213869d0*x52
      x54 = rho_a**(-0.666666666666667d0)*x53
      x55 = x0*x54
      x56 = x50 + x55
      x57 = sqrt(x56)
      x58 = 1d0/x57
      x59 = rho_a**(-0.333333333333333d0)
      x60 = mu*x59
      x61 = x37*x60
      x62 = x58*x61
      x63 = sqrt(x48 + x54)
      x64 = x49 + 0.757211d0
      x65 = sqrt(x55 + x64)
      x66 = x37*x63 - x65
      x67 = 1.10064241629821d0*x61
      x68 = 0.550321208149104d0*x60
      x69 = x37*x68
      x70 = 1d0/(x57 + x69)
      x71 = x63 + x68
      x72 = x70*x71
      x73 = log(x37*x72)
      x74 = x0*x47
      x75 = 0.151426716069345d0*x74
      x76 = x75 + 1.514422d0
      x77 = x65 + x69
      x78 = x70*x77
      x79 = log(x78)
      x80 = gamma_aa*x3
      x81 = 0.0189283395086681d0*x80 + 1.0d0
      x82 = 1d0/x81
      x83 = x80*x82
      x84 = 0.0168251906743717d0*x74 + 0.00124631042032383d0*x83 + 0.052
     &7327418253926d0
      x85 = x50**2
      x86 = 1d0/x85
      x87 = 0.0833333333333333d0*x32
      x88 = 1d0/rho_a
      x89 = x56**(-1.5d0)
      x90 = mu**3*x88*x89
      x91 = -0.825481812223657d0*x62 + x87*x90 + 1.0d0
      x92 = x86*x91
      x93 = x50**3
      x94 = x50**3.5d0
      x95 = 0.0151426716069345d0*x74 + 0.00112167937829144d0*x83 + 0.047
     &4594676428533d0
      x96 = sqrt(x46)
      x97 = x36*x37
      x98 = x96*x97
      x99 = sqrt(x64)
      x100 = 0.275160604074552d0*x98 - x99
      x101 = 2.13333333333333d0*x94
      x102 = x100*x101 - 0.888888888888889d0*x50*x95 - 0.025212152177451
     &8d0*x85 + 0.807691733333333d0*x93 + 1.26041162731059d0*x94
      x103 = 1d0/x93
      x104 = rho_a**(-1.66666666666667d0)
      x105 = 0.0189283395086681d0*x104
      x106 = 0.208333333333333d0*x32
      x107 = mu**5*x105*x22*x56**(-2.5d0) - x106*x90 + 1.03185226527957d
     &0*x62 - 1.0d0
      x108 = x103*x107
      x109 = -x102*x108 - x51*(0.0260152163368769d0*x62 - 0.047272785332
     &7222d0) + x66 *x67 + x73*x75 - x76*x79 + x84*x92 + 0.757211d0
      x110 = 0.682784063255296d0
      x111 = 2.72568088924821d0*x38
      x112 = x110*x111
      x113 = rho_a**(-3.66666666666667d0)
      x114 = x0*x113
      x115 = gamma_aa*x114
      x116 = 0.00322924538597961d0*x115
      x117 = rho_a**(-10.3333333333333d0)*x14
      x118 = 4.70452498850791d-5*x117
      x119 = rho_a**(-9.0d0)*x19
      x120 = 0.000337244097222222d0*x119
      x121 = rho_a**(-7.66666666666667d0)*x24
      x122 = 0.00160517576905591d0*x121
      x123 = rho_a**(-6.33333333333333d0)
      x124 = x123*x29
      x125 = 0.00490299508380309d0*x124
      x126 = rho_a**(-5.0d0)*x34
      x127 = 0.00710829166666667d0*x126
      x128 = x2*(x116 + x118 - x120 + x122 - x125 + x127)
      x129 = 0.151426716069345d0*x128
      x130 = x0*x42
      x131 = gamma_aa*x130
      x132 = 3.63324652777778d-6*rho_a**(-13.0d0)*x7 - 2.70095837470106d
     &-5*rho_a**( -11.6666666666667d0)*x10 + 1.95664504272581d0*rho_a**(
     & -2.33333333333333d0)*x97 - 2.51933675081772d0*x115 + 0.0003428177
     &87675372d0*x117 - 0.002119375d0*x119 + 0.0180311362159573d0*x121 -
     & 0.156234050855499d0*x124 + 0.9249d0* x126
      x133 = x44/x41**2
      x134 = x132*x133
      x135 = 0.151426716069345d0*x73
      x136 = x113*x44
      x137 = x2*(-x116 - x118 + x120 - x122 + x125 - x127)
      x138 = x132*x46
      x139 = x86*(x58*x69 - 1.0d0)
      x140 = -0.40380457618492d0*x136 + 0.151426716069345d0*x138
      x141 = x131*(0.151426716069345d0*x137 + x140)
      x142 = x104*x52
      x143 = 0.10095114404623d0*x142
      x144 = gamma_aa*x45
      x145 = x113*x144
      x146 = 0.10095114404623d0*x145
      x147 = gamma_aa*x42
      x148 = x137*x147
      x149 = 0.0378566790173362d0*x148
      x150 = gamma_aa*x134*x2
      x151 = 0.0378566790173362d0*x150
      x152 = -x143 - x146 + x149 + x151
      x153 = 1d0/x63
      x154 = -x153*x37
      x155 = 1d0/x65
      x156 = x0*x155
      x157 = x154 + x156
      x158 = 0.183440402716368d0*x40
      x159 = 0.0378566790173362d0*x128
      x160 = x143 + x146 - x151
      x161 = 1d0/x56
      x162 = x161*x32*x59
      x163 = x162*(x147*x159 + x160)
      x164 = mu*x58
      x165 = x114*x82
      x166 = x81**(-2)
      x167 = gamma_aa*x27
      x168 = x123*x166*x167
      x169 = x114*x45
      x170 = x128*x130
      x171 = 0.0168251906743717d0*x3
      x172 = x103*x84*x91
      x173 = 1d0/x71
      x174 = mu*x158
      x175 = -x149 + x160
      x176 = x0*x175*x58 + x174
      x177 = x161*x52
      x178 = rho_a**(-2.0d0)*x177
      x179 = x104*x53 + 0.30285343213869d0*x145 - 0.113570037052009d0*x1
     &50
      x180 = x56**(-2)
      x181 = x180*x22*x52*x88
      x182 = x84*x86
      x183 = x76/x77
      x184 = -0.60570686427738d0*x136 + 0.227140074104017d0*x138
      x185 = x102*x107/x50**4
      x186 = mu**4
      x187 = x180*x186
      x188 = x128*x147
      x189 = x102*x103
      x190 = -0.264996753121354d0*x128 - 0.70665800832361d0*x136 + 0.264
     &996753121354d0*x138
      x191 = x50**2.5d0
      x192 = 0.38521965630105d0*x191
      x193 = 0.0673007626974867d0*x74 + 0.541911111111111d0
      x194 = 0.0151426716069345d0*x3
      x195 = x100*x191
      x196 = 0.137580302037276d0*x98/x44
      x197 = 1d0/x99
      x198 = x167*x26
      x199 = x11*x12*x23
      x200 = x16*x17*x28
      x201 = x21*x22*x33
      x202 = x31*x32*x36
      x203 = -0.00183862315642616d0*x198 + 1.76419687069047d-5*x199 - 0.
     &000126466536458333d0*x200 + 0.000601940913395968d0*x201 + 0.002665
     &609375d0*x202 + x43
      x204 = gamma_aa*x203
      x205 = x204*x42
      x206 = x0*x46
      x207 = 0.0585877690708122d0*x198
      x208 = 1.36246744791667d-6*x13*x5*x6
      x209 = 1.0128593905129d-5*x18*x8*x9
      x210 = 0.000128556670378264d0*x199
      x211 = 0.000794765625d0*x200
      x212 = 0.00676167608098397d0*x201
      x213 = 0.3468375d0*x202
      x214 = 1d0/x36
      x215 = 0.733741891022181d0*x214*x40
      x216 = -x207 + x208 - x209 + x210 - x211 + x212 + x213 + x215 - x4
     &
      x217 = x133*x216*x80
      x218 = rho_a**(-3.0d0)
      x219 = x218*x32
      x220 = x216*x45
      x221 = x204 + x44
      x222 = -gamma_aa*x220 + x221
      x223 = x222*x42
      x224 = mu*x223
      x225 = x224*x89
      x226 = x223*x3
      x227 = 0.151426716069345d0*x226
      x228 = x3*x82
      x229 = x166*x198
      x230 = x113*x177*x22
      x231 = 0.0378566790173362d0*x226
      Ea = -x109*x112
      d1Ea(1) = -x110*(3.63424118566428d0*rho_a**0.333333333333333d0*x10
     &9 + x111*( gamma_aa*x92*(x134*x171 - 0.00332349445419687d0*x165 + 
     &6.2908231384214d-5*x168 - 0.0448671751316578d0*x169 - 0.0168251906
     &743717d0*x170) - 0.366880805432736d0*mu*x40*x66 - x108*(-gamma_aa*
     &x193*(x134*x194 - 0.00299114500877718d0*x165 + 5.66174082457926d-5
     &*x168 - 0.040380457618492d0*x169 - 0.0151426716069345d0*x170) + x1
     &01*(x1*x196*(-x128 - 2.66666666666667d0*x136 + x138) - x131*x197*(
     &-0.10095114404623d0* x136 + 0.0378566790173362d0*x138 - x159)) + 2
     &.13333333333333d0* x131*x190*x195 - 0.0252121521774518d0*x131*x50*
     &(-x129 + x140) + 0.807691733333333d0*x131*x85*(-0.227140074104017d
     &0*x128 + x184) - x131*x95*(-0.179468700526631d0*x136 + 0.067300762
     &6974867d0*x137 + 0.0673007626974867d0*x138) + 0.711111111111111d0*
     &x147*x190*x192) - 0.40380457618492d0*x115*x45*x73 - x129*x131*x73 
     &+ 0.0472727853327222d0*x131*x139*(-0.20190228809246d0*x136 + 0.075
     &7133580346725d0*x137 + 0.0757133580346725d0*x138) + x131* x185*(0.
     &227140074104017d0*x137 + x184) + x134*x135*x80 - x141* x172 - x141
     &*x79 - x152*x157*x67 - x164*x182*(0.825481812223657d0* x162*x175 +
     & x178*x87 - 0.0833333333333333d0*x181*( -0.113570037052009d0*x148 
     &+ x179) - 0.275160604074552d0*x40) + x164*x189*(-x105*x12*x186*(0.
     &50475572023115d0*x142 + 0.50475572023115d0*x145 - 0.18928339508668
     &1d0*x150 + 0.189283395086681d0*x188)/x56**3 - x106*x178 - 1.031852
     &26527957d0 *x163 + 0.208333333333333d0*x181*(x179 + 0.113570037052
     &009d0*x188 ) + 0.0315472325144469d0*x187*x2*x22 + 0.34395075509319
     &d0*x40) + 0.0472727853327222d0*x164*x51*(x158 - 0.550321208149104d
     &0*x163) + x173*x75*(-0.183440402716368d0*mu*x39 + x152*x153 + x176
     &*x72) - x183*(x152*x156 - x174 + x176*x78)))
      d1Ea(2) = -x112*(-0.0057325125848865d0*gamma_aa*x133*x173*x222*x26
     &*x37*(x154 + x32 *x58*x72) - mu*x182*x42*x89*(-0.03125d0*x219 + 0.
     &00946416975433406d0*x230)*(x144*(x207 - x208 + x209 - x210 + x211 
     &- x212 - x213 - x215 + x4) + x221) - x108*(x101*(x196*(x203 - x220
     &) - x197*x231 + 0.137580302037276d0*x214*x37*x96) + 0.188442135552
     &963d0*x192*x2*x223 - x193*(x194*x205 + 0.0151426716069345d0*x206 -
     & 0.0151426716069345d0*x217 + 0.00112167937829144d0*x228 - 2.123152
     &80921722d-5*x229) + 0.565326406658888d0*x195*x226 - 0.003817793409
     &27212d0*x226*x50 + 0.183459160162536d0*x226*x85 - 0.06730076269748
     &67d0*x226*x95) + x135*x205*x3 + x135*x206 - x135*x217 + 0.00357918
     &132119261d0*x139 *x226 - 0.0416666666666667d0*x157*x218*x224*x37 -
     & x172*x227 - x183*x231*(x155 - x58*x78) + 0.227140074104017d0*x185
     &*x226 + x189 *x225*(0.00358282036555407d0*rho_a**(-4.3333333333333
     &3d0)*x12* x187 + 0.0390625d0*x219 - 0.0236604243858352d0*x230) + 0
     &.000984849694431712d0*x219*x225*x51 - x227*x79 + x92*(x171*x205 + 
     &0.0168251906743717d0*x206 - 0.0168251906743717d0*x217 + 0.00124631
     &042032383d0*x228 - 2.35905867690803d-5*x229))
      end subroutine


C*****************************************************************************
      pure subroutine D2ESRX_wPBE(rho_a, gamma_aa, mu, Ea, d1Ea, d2Ea)
C*****************************************************************************
C   Implemented by E.R. Kjellgren.
C
C   Subroutine generated using Sympy 1.3
C   Generated: March 22, 2019
C*****************************************************************************
      implicit none
      real*8, intent(in) :: rho_a, gamma_aa, mu
      real*8, intent(out) :: Ea, d1Ea(4), d2Ea(10)
      real*8 :: x0, x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x
     &13, x14, x15, x16, x17, x18, x19, x20, x21, x22, x23, x24, x25, x2
     &6, x27, x28, x29, x30, x31, x32, x33, x34, x35, x36, x37, x38, x39
     &, x40, x41, x42, x43, x44, x45, x46, x47, x48, x49, x50, x51, x52,
     & x53, x54, x55, x56, x57, x58, x59, x60, x61, x62, x63, x64, x65, 
     &x66, x67, x68, x69, x70, x71, x72, x73, x74, x75, x76, x77, x78, x
     &79, x80, x81, x82, x83, x84, x85, x86, x87, x88, x89, x90, x91, x9
     &2, x93, x94, x95, x96, x97, x98, x99, x100, x101, x102, x103, x104
     &, x105, x106, x107, x108, x109, x110, x111, x112, x113, x114, x115
     &, x116, x117, x118, x119, x120, x121, x122, x123, x124, x125, x126
     &, x127, x128, x129, x130, x131, x132, x133, x134, x135, x136, x137
     &, x138, x139, x140, x141, x142, x143, x144, x145, x146, x147, x148
     &, x149, x150, x151, x152, x153, x154, x155, x156, x157, x158, x159
     &, x160, x161, x162, x163, x164, x165, x166, x167, x168, x169, x170
     &, x171, x172, x173, x174, x175, x176, x177, x178, x179, x180, x181
     &, x182, x183, x184, x185, x186, x187, x188, x189, x190, x191, x192
     &, x193, x194, x195, x196, x197, x198, x199, x200, x201, x202, x203
     &, x204, x205, x206, x207, x208, x209, x210, x211, x212, x213, x214
     &, x215, x216, x217, x218, x219, x220, x221, x222, x223, x224, x225
     &, x226, x227, x228, x229, x230, x231, x232, x233, x234, x235, x236
     &, x237, x238, x239, x240, x241, x242, x243, x244, x245, x246, x247
     &, x248, x249, x250, x251, x252, x253, x254, x255, x256, x257, x258
     &, x259, x260, x261, x262, x263, x264, x265, x266, x267, x268, x269
     &, x270, x271, x272, x273, x274, x275, x276, x277, x278, x279, x280
     &, x281, x282, x283, x284, x285, x286, x287, x288, x289, x290, x291
     &, x292, x293, x294, x295, x296, x297, x298, x299, x300, x301, x302
     &, x303, x304, x305, x306, x307, x308, x309, x310, x311, x312, x313
     &, x314, x315, x316, x317, x318, x319, x320, x321, x322, x323, x324
     &, x325, x326, x327, x328, x329, x330, x331, x332, x333, x334, x335
     &, x336, x337, x338, x339, x340, x341, x342, x343, x344, x345, x346
     &, x347, x348, x349, x350, x351, x352, x353, x354, x355, x356, x357
     &, x358, x359, x360, x361, x362, x363, x364, x365, x366, x367, x368
     &, x369, x370, x371, x372, x373, x374, x375, x376, x377, x378, x379
     &, x380, x381, x382, x383, x384, x385, x386, x387, x388, x389, x390
     &, x391, x392, x393, x394, x395, x396, x397, x398, x399, x400, x401
     &, x402, x403, x404, x405, x406, x407, x408, x409, x410, x411, x412
     &, x413, x414, x415, x416, x417, x418, x419, x420, x421, x422, x423
     &, x424, x425, x426, x427, x428, x429, x430, x431, x432, x433, x434
     &, x435, x436, x437, x438, x439, x440, x441, x442, x443, x444, x445
     &, x446, x447, x448, x449, x450, x451, x452, x453, x454, x455, x456
     &, x457, x458, x459, x460, x461, x462, x463, x464, x465, x466, x467
     &, x468, x469, x470, x471, x472, x473, x474, x475, x476, x477, x478
     &, x479, x480, x481, x482, x483, x484, x485, x486, x487, x488, x489
     &, x490, x491, x492, x493, x494, x495, x496, x497, x498, x499, x500
     &, x501, x502, x503, x504, x505, x506, x507, x508, x509, x510, x511
     &, x512, x513, x514, x515, x516, x517, x518, x519, x520, x521, x522
     &, x523, x524, x525, x526, x527, x528, x529, x530, x531, x532, x533
     &, x534, x535, x536, x537, x538, x539, x540, x541, x542, x543, x544
     &, x545, x546, x547, x548, x549, x550, x551, x552, x553, x554, x555
     &, x556, x557, x558, x559, x560, x561, x562, x563, x564, x565, x566
     &, x567, x568, x569, x570, x571, x572, x573, x574, x575, x576, x577
     &, x578, x579, x580, x581, x582, x583, x584, x585, x586, x587, x588
     &, x589, x590, x591, x592, x593, x594, x595, x596, x597, x598, x599
     &, x600, x601, x602, x603, x604, x605, x606, x607, x608, x609, x610
     &, x611, x612, x613, x614, x615, x616, x617, x618, x619, x620, x621
     &, x622, x623, x624, x625, x626, x627, x628, x629, x630, x631, x632
     &, x633, x634, x635, x636, x637, x638, x639, x640, x641, x642, x643
     &, x644, x645, x646, x647, x648, x649, x650, x651, x652, x653, x654
     &, x655, x656, x657, x658, x659, x660, x661, x662, x663, x664, x665
     &, x666, x667, x668, x669, x670, x671, x672, x673, x674, x675, x676
     &, x677, x678, x679, x680, x681, x682, x683, x684, x685, x686, x687
     &, x688, x689, x690, x691, x692, x693, x694, x695, x696, x697, x698
     &, x699, x700, x701, x702, x703, x704, x705, x706, x707, x708, x709
     &, x710, x711, x712, x713, x714, x715, x716, x717, x718, x719, x720
     &, x721, x722, x723, x724, x725, x726, x727, x728, x729, x730, x731
     &, x732, x733, x734, x735, x736, x737, x738, x739, x740, x741, x742
     &, x743, x744, x745, x746, x747, x748, x749, x750, x751, x752, x753
     &, x754, x755, x756, x757, x758, x759, x760, x761, x762, x763, x764
     &, x765, x766, x767, x768, x769, x770, x771, x772, x773, x774, x775
     &, x776, x777, x778, x779, x780, x781, x782, x783, x784, x785, x786
     &, x787, x788, x789, x790, x791, x792, x793, x794, x795, x796, x797
     &, x798, x799, x800, x801
      Ea = 0.0d0
      d1Ea(:) = 0.0d0
      d2Ea(:) = 0.0d0
      x0 = 0.217336917462899d0
      x1 = rho_a**2.66666666666667d0
      x2 = 1d0/x1
      x3 = x0*x2
      x4 = 0.944751281556643d0*x3
      x5 = rho_a**(-12.0d0)
      x6 = 0.00104016147329585d0
      x7 = gamma_aa**4.5d0*x6
      x8 = rho_a**(-10.6666666666667d0)
      x9 = 0.00223117693796192d0
      x10 = gamma_aa**4*x9
      x11 = rho_a**(-9.33333333333333d0)
      x12 = 0.00478594012208448d0
      x13 = gamma_aa**3.5d0
      x14 = x12*x13
      x15 = x11*x14
      x16 = rho_a**(-8.0d0)
      x17 = 0.0102659822546843d0
      x18 = gamma_aa**3
      x19 = x17*x18
      x20 = x16*x19
      x21 = rho_a**(-6.66666666666667d0)
      x22 = 0.0220208337265180d0
      x23 = gamma_aa**2.5d0
      x24 = x22*x23
      x25 = x21*x24
      x26 = gamma_aa**2
      x27 = 0.0472353356922751d0
      x28 = rho_a**5.33333333333333d0
      x29 = 1d0/x28
      x30 = x27*x29
      x31 = x26*x30
      x32 = rho_a**(-4.0d0)
      x33 = 0.101321183642338d0
      x34 = gamma_aa**1.5d0
      x35 = x33*x34
      x36 = x32*x35
      x37 = sqrt(gamma_aa)
      x38 = 0.466194077035412d0
      x39 = rho_a**1.33333333333333d0
      x40 = 1d0/x39
      x41 = x38*x40
      x42 = 1.46748378204436d0*x41
      x43 = -gamma_aa*x4 - 2.53214847628224d-6*x10*x8 + 3.67304772509327
     &d-5*x15 - 0.000264921875d0*x20 + 0.00270467043239359d0*x25 - 0.029
     &2938845354061d0*x31 + 0.231225d0*x36 + x37*x42 + 3.02770543981481d
     &-7*x5*x7 + 1.0d0
      x44 = 1d0/x43
      x45 = 0.00121096701974236d0*x3
      x46 = gamma_aa*x45 + 5.04056248768705d-6*x15 - 4.21555121527778d-5
     &*x20 + 0.000240776365358387d0*x25 - 0.000919311578213079d0*x31 + 0
     &.00177707291666667d0*x36
      x47 = x44*x46
      x48 = x2*x47
      x49 = 0.0757133580346725d0*x48
      x50 = gamma_aa*x49
      x51 = x0*x50
      x52 = x51 + 0.60965d0
      x53 = 1d0/x52
      x54 = mu**2
      x55 = 0.30285343213869d0*x54
      x56 = rho_a**(-0.666666666666667d0)*x55
      x57 = x0*x56
      x58 = x52 + x57
      x59 = sqrt(x58)
      x60 = 1d0/x59
      x61 = mu*x60
      x62 = rho_a**(-0.333333333333333d0)
      x63 = x38*x62
      x64 = x61*x63
      x65 = x50 + x56
      x66 = sqrt(x65)
      x67 = x51 + 0.757211d0
      x68 = x57 + x67
      x69 = sqrt(x68)
      x70 = x38*x66 - x69
      x71 = mu*x70
      x72 = 1.10064241629821d0*x63
      x73 = 0.550321208149104d0*x62
      x74 = mu*x73
      x75 = x38*x74
      x76 = x59 + x75
      x77 = 1d0/x76
      x78 = x66 + x74
      x79 = x77*x78
      x80 = x38*x79
      x81 = log(x80)
      x82 = gamma_aa*x0
      x83 = x48*x82
      x84 = 0.151426716069345d0*x83
      x85 = x84 + 1.514422d0
      x86 = x69 + x75
      x87 = x77*x86
      x88 = log(x87)
      x89 = gamma_aa*x3
      x90 = 0.0189283395086681d0*x89 + 1.0d0
      x91 = 1d0/x90
      x92 = x89*x91
      x93 = 0.0168251906743717d0*x83 + 0.00124631042032383d0*x92 + 0.052
     &7327418253926d0
      x94 = x52**2
      x95 = 1d0/x94
      x96 = 0.0833333333333333d0*x33
      x97 = x58**(-1.5d0)
      x98 = 1d0/rho_a
      x99 = mu**3*x97*x98
      x100 = -0.825481812223657d0*x64 + x96*x99 + 1.0d0
      x101 = x100*x95
      x102 = x52**3
      x103 = x52**3.5d0
      x104 = 0.0151426716069345d0*x83 + 0.00112167937829144d0*x92 + 0.04
     &74594676428533d0
      x105 = x38*sqrt(x48)
      x106 = x105*x37
      x107 = sqrt(x67)
      x108 = 0.275160604074552d0*x106 - x107
      x109 = 2.13333333333333d0*x103
      x110 = 0.807691733333333d0*x102 + 1.26041162731059d0*x103 - 0.8888
     &88888888889d0 *x104*x52 + x108*x109 - 0.0252121521774518d0*x94
      x111 = 1d0/x102
      x112 = rho_a**1.66666666666667d0
      x113 = 1d0/x112
      x114 = 0.0189283395086681d0*x113
      x115 = x22*x58**(-2.5d0)
      x116 = 0.208333333333333d0*x33
      x117 = mu**5*x114*x115 - x116*x99 + 1.03185226527957d0*x64 - 1.0d0
     &
      x118 = x111*x117
      x119 = x101*x93 - x110*x118 - x53*(0.0260152163368769d0*x64 - 0.04
     &72727853327222d0) + x71*x72 + x81*x84 - x85*x88 + 0.757211d0
      x120 = 0.682784063255296d0
      x121 = 2.72568088924821d0*x39
      x122 = x120*x121
      x123 = rho_a**0.333333333333333d0
      x124 = 3.63424118566428d0*x123
      x125 = mu*x41
      x126 = 0.366880805432736d0*x125
      x127 = 0.40380457618492d0*x81
      x128 = rho_a**(-3.66666666666667d0)
      x129 = x0*x128
      x130 = gamma_aa*x129
      x131 = x130*x47
      x132 = x0*x44
      x133 = gamma_aa*x132
      x134 = 0.00322924538597961d0*x129
      x135 = gamma_aa*x134
      x136 = rho_a**(-10.3333333333333d0)
      x137 = x136*x14
      x138 = 4.70452498850791d-5*x137
      x139 = rho_a**(-9.0d0)
      x140 = x139*x19
      x141 = 0.000337244097222222d0*x140
      x142 = rho_a**(-7.66666666666667d0)
      x143 = x142*x24
      x144 = 0.00160517576905591d0*x143
      x145 = rho_a**(-6.33333333333333d0)
      x146 = x26*x27
      x147 = x145*x146
      x148 = 0.00490299508380309d0*x147
      x149 = rho_a**(-5.0d0)
      x150 = x149*x35
      x151 = 0.00710829166666667d0*x150
      x152 = x135 + x138 - x141 + x144 - x148 + x151
      x153 = x152*x2
      x154 = 0.151426716069345d0*x153
      x155 = x154*x81
      x156 = 2.51933675081772d0*x129
      x157 = rho_a**(-13.0d0)
      x158 = x157*x7
      x159 = rho_a**(-11.6666666666667d0)
      x160 = x10*x159
      x161 = rho_a**(-2.33333333333333d0)*x38
      x162 = x161*x37
      x163 = -gamma_aa*x156 + 0.000342817787675372d0*x137 - 0.002119375d
     &0*x140 + 0.0180311362159573d0*x143 - 0.156234050855499d0*x147 + 0.
     &9249d0* x150 + 3.63324652777778d-6*x158 - 2.70095837470106d-5*x160
     & + 1.95664504272581d0*x162
      x164 = x43**(-2)
      x165 = x164*x46
      x166 = x163*x165
      x167 = 0.151426716069345d0*x81
      x168 = x167*x89
      x169 = x128*x46
      x170 = 0.20190228809246d0*x169
      x171 = -x135 - x138 + x141 - x144 + x148 - x151
      x172 = x171*x2
      x173 = 0.0757133580346725d0*x172
      x174 = x163*x48
      x175 = -x170 + x173 + 0.0757133580346725d0*x174
      x176 = x60*x75 - 1.0d0
      x177 = x176*x95
      x178 = 0.0472727853327222d0*x177
      x179 = x133*x178
      x180 = 0.151426716069345d0*x172
      x181 = 0.40380457618492d0*x169
      x182 = 0.151426716069345d0*x174 - x181
      x183 = x180 + x182
      x184 = x132*x88
      x185 = gamma_aa*x184
      x186 = mu*x72
      x187 = x113*x54
      x188 = 0.10095114404623d0*x187
      x189 = x128*x47
      x190 = gamma_aa*x189
      x191 = 0.10095114404623d0*x190
      x192 = 0.0378566790173362d0*x172
      x193 = gamma_aa*x44
      x194 = x192*x193
      x195 = x163*x2
      x196 = x165*x195
      x197 = gamma_aa*x196
      x198 = 0.0378566790173362d0*x197
      x199 = -x188 - x191 + x194 + x198
      x200 = 1d0/x66
      x201 = x200*x38
      x202 = -x201
      x203 = 1d0/x69
      x204 = x0*x203
      x205 = x202 + x204
      x206 = x199*x205
      x207 = 0.183440402716368d0*x41
      x208 = 0.0378566790173362d0*x153
      x209 = x193*x208
      x210 = x188 + x191 - x198
      x211 = x209 + x210
      x212 = 1d0/x58
      x213 = x212*x33
      x214 = x211*x213
      x215 = x207 - x214*x73
      x216 = 0.0472727853327222d0*x53
      x217 = x216*x61
      x218 = x90**(-2)
      x219 = gamma_aa*x27
      x220 = x145*x219
      x221 = x218*x220
      x222 = x129*x91
      x223 = x129*x47
      x224 = x153*x44
      x225 = x0*x224
      x226 = 0.0168251906743717d0*x3
      x227 = x166*x226 - 0.00332349445419687d0*x222 - 0.0448671751316578
     &d0*x223 - 0.0168251906743717d0*x225
      x228 = 6.2908231384214d-5*x221 + x227
      x229 = gamma_aa*x101
      x230 = x100*x111
      x231 = x230*x93
      x232 = x133*x183
      x233 = 1d0/x78
      x234 = 0.183440402716368d0*mu*x40
      x235 = x199*x200
      x236 = mu*x207
      x237 = -x194 + x210
      x238 = x0*x60
      x239 = x237*x238
      x240 = x236 + x239
      x241 = -x234 + x235 + x240*x79
      x242 = x233*x241
      x243 = x242*x84
      x244 = rho_a**(-2.0d0)
      x245 = x212*x54
      x246 = x244*x245
      x247 = x213*x237
      x248 = 0.825481812223657d0*x62
      x249 = 0.0833333333333333d0*x98
      x250 = 0.113570037052009d0*x172
      x251 = x193*x250
      x252 = 0.113570037052009d0*gamma_aa
      x253 = x196*x252
      x254 = x113*x55 + 0.30285343213869d0*x190 - x253
      x255 = -x251 + x254
      x256 = x58**(-2)
      x257 = x22*x256
      x258 = x255*x257
      x259 = x258*x54
      x260 = x246*x96 + x247*x248 - x249*x259 - 0.275160604074552d0*x41
      x261 = x93*x95
      x262 = x261*x61
      x263 = x199*x204
      x264 = -x236 + x240*x87 + x263
      x265 = 1d0/x86
      x266 = x265*x85
      x267 = x264*x266
      x268 = 0.227140074104017d0*x172
      x269 = -0.60570686427738d0*x169 + 0.227140074104017d0*x174
      x270 = x268 + x269
      x271 = x52**(-4)
      x272 = x117*x271
      x273 = x110*x272
      x274 = x133*x273
      x275 = mu**4
      x276 = x256*x275
      x277 = x22*x276
      x278 = 1.03185226527957d0*x62
      x279 = x114*x275
      x280 = 0.189283395086681d0*x153
      x281 = 0.189283395086681d0*x197
      x282 = 0.50475572023115d0*x187 + 0.50475572023115d0*x190 - x281
      x283 = x58**(-3)
      x284 = x12*x283
      x285 = x284*(x193*x280 + x282)
      x286 = x224*x252 + x254
      x287 = x257*x286
      x288 = 0.208333333333333d0*x98
      x289 = x288*x54
      x290 = -x116*x246 + 0.0315472325144469d0*x2*x277 - x214*x278 - x27
     &9*x285 + x287 *x289 + 0.34395075509319d0*x41
      x291 = x110*x111
      x292 = x291*x61
      x293 = -x154 + x182
      x294 = 0.227140074104017d0*x153
      x295 = x269 - x294
      x296 = 0.807691733333333d0*x94
      x297 = x132*x296
      x298 = gamma_aa*x297
      x299 = x52**2.5d0
      x300 = 0.264996753121354d0*x153
      x301 = -0.70665800832361d0*x169 + 0.264996753121354d0*x174 - x300
      x302 = x299*x301
      x303 = 0.385219656301050d0
      x304 = 0.711111111111111d0*x303
      x305 = x193*x304
      x306 = 0.0673007626974867d0*x172
      x307 = -0.179468700526631d0*x169 + 0.0673007626974867d0*x174 + x30
     &6
      x308 = x104*x132
      x309 = gamma_aa*x308
      x310 = 0.0151426716069345d0*x3
      x311 = x166*x310 - 0.00299114500877718d0*x222 - 0.040380457618492d
     &0*x223 - 0.0151426716069345d0*x225
      x312 = 5.66174082457926d-5*x221 + x311
      x313 = 0.0673007626974867d0*gamma_aa
      x314 = x0*x48
      x315 = x313*x314 + 0.541911111111111d0
      x316 = gamma_aa*x315
      x317 = 2.13333333333333d0*x108
      x318 = x133*x317
      x319 = -x153 - 2.66666666666667d0*x169 + x174
      x320 = 1d0/x46
      x321 = 0.137580302037276d0*x106
      x322 = x320*x321
      x323 = x1*x322
      x324 = x319*x323
      x325 = 0.10095114404623d0*x169
      x326 = 0.0378566790173362d0*x174
      x327 = -x325 + x326
      x328 = -x208 + x327
      x329 = x44/x107
      x330 = x329*x82
      x331 = x324 - x328*x330
      x332 = x109*x331 - 0.0252121521774518d0*x133*x293*x52 + x295*x298 
     &+ x302*x305 + x302*x318 - x307*x309 - x312*x316
      x333 = -x118*x332 - x126*x70 - x127*x131 - x133*x155 + x166*x168 +
     & x175*x179 - x183*x185 - x186*x206 + x215*x217 + x228*x229 - x231*
     &x232 + x243 - x260*x262 - x267 + x270*x274 + x290*x292
      x334 = 0.0585877690708122d0*x30
      x335 = gamma_aa*x334
      x336 = x13*x6
      x337 = x336*x5
      x338 = 1.36246744791667d-6*x337
      x339 = x18*x9
      x340 = x339*x8
      x341 = 1.0128593905129d-5*x340
      x342 = x12*x23
      x343 = x11*x342
      x344 = 0.000128556670378264d0*x343
      x345 = x17*x26
      x346 = x16*x345
      x347 = 0.000794765625d0*x346
      x348 = x22*x34
      x349 = x21*x348
      x350 = 0.00676167608098397d0*x349
      x351 = x33*x37
      x352 = x32*x351
      x353 = 0.3468375d0*x352
      x354 = 1d0/x37
      x355 = 0.733741891022181d0*x354*x41
      x356 = -x335 + x338 - x341 + x344 - x347 + x350 + x353 + x355 - x4
     &
      x357 = x356*x47
      x358 = gamma_aa*x357
      x359 = 0.00183862315642616d0*x30
      x360 = -gamma_aa*x359 + 1.76419687069047d-5*x343 - 0.0001264665364
     &58333d0*x346 + 0.000601940913395968d0*x349 + 0.002665609375d0*x352
     & + x45
      x361 = gamma_aa*x360
      x362 = x361 + x46
      x363 = -x358 + x362
      x364 = x3*x44
      x365 = x104*x364
      x366 = 0.0673007626974867d0*x365
      x367 = x3*x91
      x368 = gamma_aa*x30
      x369 = x218*x368
      x370 = x361*x44
      x371 = x165*x356
      x372 = x371*x89
      x373 = x310*x370 + 0.0151426716069345d0*x314 + 0.00112167937829144
     &d0*x367 - 2.12315280921722d-5*x369 - 0.0151426716069345d0*x372
      x374 = x105*x354
      x375 = 0.137580302037276d0*x374
      x376 = -x357 + x360
      x377 = x320*x376
      x378 = x321*x377
      x379 = x3*x329
      x380 = 0.0378566790173362d0*x363
      x381 = x375 + x378 - x379*x380
      x382 = x2*x299
      x383 = x363*x44
      x384 = x363*x364
      x385 = x299*x384
      x386 = 0.565326406658888d0*x385
      x387 = x108*x386 + x109*x381 + 0.188442135552963d0*x303*x382*x383 
     &- x315*x373 - 0.00381779340927212d0*x384*x52 + 0.183459160162536d0
     &*x384*x94
      x388 = x364*x88
      x389 = 0.151426716069345d0*x388
      x390 = 0.00357918132119261d0*x177
      x391 = rho_a**(-3.0d0)
      x392 = mu*x205*x38
      x393 = 0.0416666666666667d0*x392
      x394 = x33*x391
      x395 = rho_a**(-4.33333333333333d0)
      x396 = x12*x276
      x397 = x395*x396
      x398 = x22*x245
      x399 = x128*x398
      x400 = 0.0390625d0*x394 + 0.00358282036555407d0*x397 - 0.023660424
     &3858352d0* x399
      x401 = mu*x97
      x402 = x383*x401
      x403 = 0.03125d0*x394
      x404 = 0.00946416975433406d0*x399
      x405 = x335 - x338 + x341 - x344 + x347 - x350 - x353 - x355 + x4
      x406 = x405*x47
      x407 = gamma_aa*x406
      x408 = x362 + x407
      x409 = x401*x44
      x410 = x261*x409
      x411 = x408*x410
      x412 = 0.227140074104017d0*x273
      x413 = 0.000984849694431712d0*x394*x53
      x414 = x33*x60
      x415 = x414*x79
      x416 = x202 + x415
      x417 = x233*x38
      x418 = x165*x417
      x419 = x416*x418
      x420 = 0.0057325125848865d0*gamma_aa
      x421 = x226*x370 + 0.0168251906743717d0*x314 + 0.00124631042032383
     &d0*x367 - 2.35905867690803d-5*x369 - 0.0168251906743717d0*x372
      x422 = x167*x3
      x423 = x60*x87
      x424 = x203 - x423
      x425 = x266*x44
      x426 = x101*x421 + x167*x314 - x167*x372 - x3*x380*x424*x425 + x37
     &0*x422
      x427 = rho_a**(-4.66666666666667d0)
      x428 = x0*x427
      x429 = gamma_aa*x428
      x430 = x128*x152
      x431 = 0.80760915236984d0*x430
      x432 = x133*x81
      x433 = rho_a**(-11.3333333333333d0)*x14
      x434 = rho_a**(-10.0d0)
      x435 = x19*x434
      x436 = rho_a**(-8.66666666666667d0)*x24
      x437 = rho_a**(-7.33333333333333d0)
      x438 = x146*x437
      x439 = rho_a**(-6.0d0)*x35
      x440 = x2*(0.0118405664152586d0*x429 + 0.000486134248812484d0*x433
     & - 0.003035196875d0*x435 + 0.012306347562762d0*x436 - 0.0310523021
     &974196d0*x438 + 0.0355414583333333d0*x439)
      x441 = 0.151426716069345d0*x440
      x442 = 0.80760915236984d0*x169
      x443 = x163*x164
      x444 = x443*x82
      x445 = 0.30285343213869d0*x81
      x446 = 9.23756808633162d0*x429
      x447 = 4.72322048611111d-5*rho_a**(-14.0d0)*x7
      x448 = 0.00031511181038179d0*rho_a**(-12.6666666666667d0)*x10
      x449 = 0.00354245047264551d0*x433
      x450 = 0.019074375d0*x435
      x451 = 0.138238710989006d0*x436
      x452 = 0.989482322084827d0*x438
      x453 = 4.6245d0*x439
      x454 = 4.56550509969357d0*rho_a**(-3.33333333333333d0)*x37*x38
      x455 = -x446 + x447 - x448 + x449 - x450 + x451 - x452 + x453 + x4
     &54
      x456 = x165*x455
      x457 = -5.03867350163543d0*x130 + 0.000685635575350743d0*x137 - 0.
     &00423875d0* x140 + 0.0360622724319145d0*x143 - 0.312468101710998d0
     &*x147 + 1.8498d0*x150 + 7.26649305555556d-6*x158 - 5.4019167494021
     &1d-5* x160 + 3.91329008545163d0*x162
      x458 = x46/x43**3
      x459 = x457*x458
      x460 = x163*x459
      x461 = x146*x164
      x462 = x183*x461
      x463 = x111*x176
      x464 = x175*x463
      x465 = x270*x271
      x466 = x100*x93
      x467 = mu*x238
      x468 = x193*x467
      x469 = x215*x95
      x470 = 2.0d0*x132*x26
      x471 = x228*x230
      x472 = x427*x46
      x473 = x128*x171
      x474 = x163*x189
      x475 = x163*x44
      x476 = x446 - x447 + x448 - x449 + x450 - x451 + x452 - x453 - x45
     &4
      x477 = x196*x457
      x478 = x476*x48
      x479 = 0.30285343213869d0*x172
      x480 = x441 + 1.48061677934471d0*x472 - 0.80760915236984d0*x474 + 
     &0.151426716069345d0*x477
      x481 = -0.80760915236984d0*x473 + x475*x479 + 0.151426716069345d0*
     &x478 + x480
      x482 = x428*x91
      x483 = x218*x219*x437
      x484 = x90**(-3)
      x485 = x345*x434*x484
      x486 = x428*x47
      x487 = x132*x430
      x488 = x132*x440
      x489 = x0*x443
      x490 = x169*x489
      x491 = x153*x489
      x492 = 0.30285343213869d0*x242
      x493 = x264*x265
      x494 = x260*x61
      x495 = x228*x95
      x496 = 2.0d0*x468
      x497 = x111*x93
      x498 = x260*x497
      x499 = x110*x117/x52**5
      x500 = x270*x499
      x501 = x78**(-2)
      x502 = x241*x501
      x503 = x85/x86**2
      x504 = x264*x503
      x505 = x77*(-x199*x238 + x236)
      x506 = x172*x475
      x507 = 0.227140074104017d0*x440 + 2.22092516901706d0*x472 - 1.2114
     &1372855476d0* x474 + 0.227140074104017d0*x477
      x508 = 0.244587203621824d0*x161
      x509 = x247*x40
      x510 = x237*x258
      x511 = x2*x54
      x512 = 0.168251906743717d0*x511
      x513 = gamma_aa*x47
      x514 = x427*x513
      x515 = 0.370154194836177d0*x514
      x516 = x193*x473
      x517 = 0.20190228809246d0*x516
      x518 = 0.0378566790173362d0*x440
      x519 = x193*x518
      x520 = gamma_aa*x443
      x521 = x170*x520
      x522 = x173*x520
      x523 = gamma_aa*x2
      x524 = 0.0378566790173362d0*x523
      x525 = x165*x476
      x526 = x524*x525
      x527 = x195*x459
      x528 = 0.0378566790173362d0*gamma_aa
      x529 = x527*x528
      x530 = x213*(-x512 - x515 + x517 - x519 + x521 - x522 - x526 - x52
     &9)
      x531 = x110*x290
      x532 = x65**(-1.5d0)
      x533 = 0.466194077035412d0
      x534 = x237*x533
      x535 = x532*x534
      x536 = x68**(-1.5d0)
      x537 = x237*x27
      x538 = x199*x537
      x539 = x512 + x515 + x519 - x521 + x529
      x540 = -x517 + x522 + x526 + x539
      x541 = x133*x270
      x542 = x272*x332
      x543 = x245*x394
      x544 = x249*x54
      x545 = 0.50475572023115d0*x511
      x546 = 1.11046258450853d0*x514
      x547 = 0.60570686427738d0*x516
      x548 = x193*x440
      x549 = 0.113570037052009d0*x548
      x550 = x169*x520
      x551 = 0.60570686427738d0*x550
      x552 = x268*x520
      x553 = x2*x252*x525
      x554 = x252*x527
      x555 = x290*x61
      x556 = x111*x332
      x557 = mu*x508
      x558 = 2.0d0*x240
      x559 = x77*(x234 - x235)
      x560 = 0.20190228809246d0*x430
      x561 = -0.0757133580346725d0*x153*x520 + x193*x560 - x456*x524 + x
     &539
      x562 = x76**(-2)
      x563 = x240*x562
      x564 = x563*(x126 + 2.0d0*x239)
      x565 = x238*x540 + x538*x97 + x557
      x566 = 0.151426716069345d0*gamma_aa
      x567 = x77*(x236 - x263)
      x568 = 0.189283395086681d0*x172
      x569 = x193*x568
      x570 = x282 - x569
      x571 = x284*x570
      x572 = 0.00104016147329585d0
      x573 = x570*x572
      x574 = x293*x461
      x575 = x295*x52
      x576 = x52**1.5d0
      x577 = 0.0837224526465878d0*x576
      x578 = 0.50475572023115d0*x169
      x579 = 0.189283395086681d0*x174
      x580 = x301*(-x280 - x578 + x579)
      x581 = x455*x48
      x582 = x163*x224
      x583 = 0.00190889670463606d0*x83 + 0.0153705885749835d0
      x584 = x132*x583
      x585 = x299*(1.41331601664722d0*x430 + 0.264996753121354d0*x440 + 
     &2.59107936385324d0*x472 - 1.41331601664722d0*x474 + 0.264996753121
     &354d0*x477 - 0.264996753121354d0*x581 - 0.529993506242707d0*x582)
      x586 = x133*x302
      x587 = x319*x320
      x588 = x46**(-2)
      x589 = 0.0687901510186381d0*x106*x588
      x590 = x321*x588
      x591 = x1*x319
      x592 = x192 + x327
      x593 = x67**(-1.5d0)
      x594 = x328*x593
      x595 = -x366*x408 + x387
      x596 = x364*x408
      x597 = x391*x408
      x598 = x291*x409
      x599 = x400*x408
      x600 = x403 - x404
      x601 = x408*x409
      x602 = 0.151426716069345d0*x596
      x603 = x29*x408
      x604 = x419*x603
      x605 = x129*x370
      x606 = 0.00980599016760618d0*x220
      x607 = x136*x342
      x608 = 0.000164658374597777d0*x607
      x609 = x139*x345
      x610 = 0.00101173229166667d0*x609
      x611 = x142*x348
      x612 = 0.00401293942263979d0*x611
      x613 = x149*x351
      x614 = 0.0106624375d0*x613
      x615 = x2*(x134 - x606 + x608 - x610 + x612 + x614)
      x616 = gamma_aa*x615
      x617 = 0.151426716069345d0*x616
      x618 = x164*x356
      x619 = x618*x82
      x620 = x361*x443
      x621 = -x156 + 1.6349609375d-5*x157*x336 - 0.000108038334988042d0*
     &x159*x339 + 0.978322521362907d0*x161*x354 - 0.312468101710998d0*x2
     &20 + 0.0011998622568638d0*x607 - 0.006358125d0*x609 + 0.0450778405
     &398931d0*x611 + 1.38735d0*x613
      x622 = x165*x89
      x623 = x621*x622
      x624 = 1.88950256311329d0*x3
      x625 = 0.117175538141624d0*x368
      x626 = 2.72493489583333d-6*x337
      x627 = 2.02571878102579d-5*x340
      x628 = 0.000257113340756529d0*x343
      x629 = 0.00158953125d0*x346
      x630 = 0.0135233521619679d0*x349
      x631 = 0.693675d0*x352
      x632 = x354*x42
      x633 = -x624 - x625 + x626 - x627 + x628 - x629 + x630 + x631 + x6
     &32
      x634 = x458*x633
      x635 = x168*x634
      x636 = x408*x44
      x637 = x22*x401*x95
      x638 = x2*x408
      x639 = x164*x219
      x640 = x638*x639
      x641 = x416*x417*x603
      x642 = x164*x641
      x643 = x458*x641
      x644 = x401*x600
      x645 = x408*x644
      x646 = x164*x82
      x647 = x271*x466
      x648 = x44*x467*x638
      x649 = x636*x89
      x650 = x230*x421
      x651 = x2*x380
      x652 = x265*x424
      x653 = 0.151426716069345d0*x242
      x654 = x2*(-x134 + x606 - x608 + x610 - x612 - x614)
      x655 = gamma_aa*x654
      x656 = x128*x361
      x657 = x195*x370
      x658 = x48*x621
      x659 = gamma_aa*x658
      x660 = x193*x405
      x661 = x128*x407
      x662 = x624 + x625 - x626 + x627 - x628 + x629 - x630 - x631 - x63
     &2
      x663 = x197*x662
      x664 = -0.40380457618492d0*x656 + 0.151426716069345d0*x657 + 0.151
     &426716069345d0*x659
      x665 = x180*x660 + x183 + 0.151426716069345d0*x655 - 0.40380457618
     &492d0*x661 + 0.151426716069345d0*x663 + x664
      x666 = x421*x95
      x667 = x200*x29
      x668 = 0.0057325125848865d0*x165*x363
      x669 = x484*x609
      x670 = x132*x616
      x671 = x169*x619
      x672 = x153*x619
      x673 = x634*x89
      x674 = x163*x673
      x675 = x2*x204
      x676 = x380*x675
      x677 = x110*x401
      x678 = x60*x77
      x679 = x2*x238
      x680 = x380*x679
      x681 = x32*x33
      x682 = x408*x681
      x683 = x212*x22
      x684 = x237*x597*x683
      x685 = 0.10095114404623d0*x656
      x686 = 0.0378566790173362d0*x655
      x687 = 0.10095114404623d0*x661
      x688 = 0.0378566790173362d0*x657
      x689 = x194*x405
      x690 = 0.0378566790173362d0*x659
      x691 = x198*x662
      x692 = x33*(-x192 + x325 - x326 + x685 - x686 + x687 - x688 - x689
     & - x690 - x691)
      x693 = -0.60570686427738d0*x656 + 0.227140074104017d0*x657 + 0.227
     &140074104017d0*x659
      x694 = 0.0378566790173362d0*x638
      x695 = x199*x694
      x696 = x532*x533
      x697 = x27*x695
      x698 = -x685 + x688 + x690
      x699 = x592 + x686 - x687 + x689 + x691 + x698
      x700 = x186*x44
      x701 = x409*x599
      x702 = x272*x595
      x703 = x398*x408*x427
      x704 = x12*x256*x54
      x705 = x128*x255*x408*x704
      x706 = x62*x692
      x707 = x398*(0.30285343213869d0*x169 - 0.113570037052009d0*x174 - 
     &x250 - x251* x405 - x252*x658 - x253*x662 - 0.113570037052009d0*x6
     &55 + 0.30285343213869d0*x656 - 0.113570037052009d0*x657 + 0.302853
     &43213869d0*x661)
      x708 = x111*x595
      x709 = x240*x77
      x710 = 0.0757133580346725d0*x363*x563
      x711 = x128*x358
      x712 = x165*x633
      x713 = x195*x712
      x714 = x209*x356 + x328 - x528*x713 - 0.0378566790173362d0*x616 + 
     &x698 + 0.10095114404623d0*x711
      x715 = x0*x699 - x212*x697
      x716 = 0.151426716069345d0*x418*x523
      x717 = x0*x714
      x718 = x275*x283
      x719 = x164*x363
      x720 = x2*x219*x719
      x721 = x108*x576
      x722 = x193*x356
      x723 = gamma_aa*x713
      x724 = x299*(x300*x722 + x301 - 0.264996753121354d0*x616 - 0.70665
     &800832361d0* x656 + 0.264996753121354d0*x657 + 0.264996753121354d0
     &*x659 + 0.70665800832361d0*x711 - 0.264996753121354d0*x723)
      x725 = x304*x44
      x726 = x313*x662
      x727 = x360*x590
      x728 = x356*x44
      x729 = 2.66666666666667d0*x128
      x730 = x360*x44
      x731 = gamma_aa*x16
      x732 = x17*x731
      x733 = x11*x12*x34
      x734 = x21*x22*x37
      x735 = x354*x681
      x736 = -x359 - 0.000252933072916667d0*x732 + 4.41049217672617d-5*x
     &733 + 0.000902911370093952d0*x734 + 0.0013328046875d0*x735
      x737 = gamma_aa*x736
      x738 = 0.151426716069345d0*x737
      x739 = x3*x730
      x740 = 0.00158953125d0*x732
      x741 = 4.76863606770833d-6*x23*x5*x6
      x742 = 3.03857817153869d-5*x26*x8*x9
      x743 = 0.000321391675945661d0*x733
      x744 = 0.010142514121476d0*x734
      x745 = 0.17341875d0*x735
      x746 = 1d0/x34
      x747 = 0.36687094551109d0*x41*x746
      x748 = -x334 - x740 + x741 - x742 + x743 + x744 + x745 - x747
      x749 = x622*x748
      x750 = x3*x361*x618
      x751 = x3*x371
      x752 = rho_a**(-5.66666666666667d0)
      x753 = x408**2
      x754 = x164*x753
      x755 = x752*x754
      x756 = x30*x754
      x757 = x3*x754
      x758 = x334 + x740 - x741 + x742 - x743 - x744 - x745 + x747
      x759 = x513*x758
      x760 = x370*x405
      x761 = x165*x405
      x762 = x662*x761
      x763 = gamma_aa*x762
      x764 = 9.16863795339332d-5*x3 + 1.33573269314236d-6*x343 - 9.57520
     &615427475d-6* x346 + 4.55749678916667d-5*x349 + 0.0002018222369899
     &55d0*x352 - 0.000139208333333333d0*x368 + 0.0378566790173362d0*x73
     &7
      x765 = 0.0757133580346725d0*x406 + 0.0378566790173362d0*x759 + 0.0
     &757133580346725d0*x760 + 0.0378566790173362d0*x763 + x764
      x766 = x394*x765
      x767 = 0.000366745518135733d0*x3 + 5.34293077256945d-6*x343 - 3.83
     &00824617099d-5*x346 + 0.000182299871566667d0*x349 + 0.000807288947
     &959818d0*x352 - 0.000556833333333333d0*x368 + x738
      x768 = 0.30285343213869d0*x406 + x566*x762 + 0.151426716069345d0*x
     &759 + 0.30285343213869d0*x760 + x767
      x769 = x363**2
      x770 = x164*x769
      x771 = x30*x770
      x772 = 0.00143312814622163d0*x424*x771
      x773 = x30*x408
      x774 = x719*x773
      x775 = 0.000217013888888889d0*x363*x408*x416*x458*x731
      x776 = x218*x30
      x777 = x484*x732
      x778 = x364*x737
      x779 = x356*x673
      x780 = 0.000550118277203599d0*x3 + 8.01439615885417d-6*x343 - 5.74
     &512369256485d-5*x346 + 0.00027344980735d0*x349 + 0.001210933421939
     &73d0*x352 - 0.00083525d0*x368 + 0.227140074104017d0*x737
      x781 = x44*x753
      x782 = 0.00143312814622163d0*x781
      x783 = x2*x765
      x784 = x29*x782
      x785 = -x201*x783 + x696*x784
      x786 = x683*x781
      x787 = x752*x786
      x788 = x145*x704*x781
      x789 = x47*x758
      x790 = x399*(x252*x762 + x252*x789 + 0.0002750591386018d0*x3 + 4.0
     &0719807942708d-6*x343 - 2.87256184628243d-5*x346 + 0.0001367249036
     &75d0*x349 + 0.000605466710969864d0*x352 - 0.000417625d0*x368 + 0.2
     &27140074104017d0*x406 + 0.113570037052009d0*x737 + 0.2271400741040
     &17d0*x760)
      x791 = 0.00286625629244325d0*x77
      x792 = 0.00286625629244325d0*x562
      x793 = 0.00143312814622163d0*x383*x773
      x794 = x30*x44*x769
      x795 = x47*x748
      x796 = gamma_aa*x795
      x797 = x356*x370
      x798 = x356*x712
      x799 = gamma_aa*x798
      x800 = -0.0757133580346725d0*x357 + x764 - 0.0378566790173362d0*x7
     &96 - 0.0757133580346725d0*x797 + 0.0378566790173362d0*x799
      x801 = 0.000641804656737533d0*x3 + 9.35012885199653d-6*x343 - 6.70
     &264430799233d-5*x346 + 0.000319024775241667d0*x349 + 0.00141275565
     &892968d0*x352 - 0.529993506242707d0*x357 - 0.000974458333333333d0*
     &x368 + 0.264996753121354d0*x737 - 0.264996753121354d0*x796 - 0.529
     &993506242707d0*x797 + 0.264996753121354d0*x799
      Ea = -x119*x122
      d1Ea(1) = -x120*(x119*x124 + x121*x333)
      d1Ea(2) = -x122*(-x118*(-x363*x366 + x387) - 0.151426716069345d0*x
     &231*x384 - x29* x363*x419*x420 + x291*x400*x402 - x363*x389 - x383
     &*x391*x393 + x384*x390 + x384*x412 + x402*x413 - x411*(-x403 + x40
     &4) + x426)
      d2Ea(1) = -x120*(1.21141372855476d0*rho_a**(-0.666666666666667d0)*
     &x119 + x121*( -2.0d0*gamma_aa*x494*x495 - x118*(-gamma_aa*x584*(x4
     &31 + x480 - 0.151426716069345d0*x581 - 0.30285343213869d0*x582) + 
     &x109*( 0.366880805432736d0*x106*x112*x587 + x152*x590*x591 + x28*x
     &319**2 *x589 + x323*(5.33333333333333d0*x430 + x440 + 9.7777777777
     &7778d0 *x472 - 5.33333333333333d0*x474 + x477 - x581 - 2.0d0*x582)
     & - x324*x475 - x330*(0.370154194836177d0*x472 - 0.20190228809246d0
     &* x474 + 0.0378566790173362d0*x477 + x518 + x560 - 0.0378566790173
     &362d0*x581 - 0.0757133580346725d0*x582) + x461* x592*x594) + 0.711
     &111111111111d0*x164*x26*x577*x580 + x298*( 1.21141372855476d0*x430
     & + x507 - 0.227140074104017d0*x581 - 0.454280148208035d0*x582) + x
     &305*x585 - x307*x312*x470 - x309*( 0.0673007626974867d0*x440 + 0.6
     &58051901930981d0*x472 - 0.358937401053262d0*x473 - 0.3589374010532
     &62d0*x474 + 0.0673007626974867d0*x477 + 0.0673007626974867d0*x478 
     &+ 0.134601525394973d0*x506) - x316*(-x310*x456 + x310*x460 + 0.010
     &9675316988497d0*x482 - 0.000509556674212134d0*x483 + 5.71559213534
     &522d-6*x485 + 0.148061677934471d0*x486 + 0.080760915236984d0*x487 
     &+ 0.0151426716069345d0*x488 - 0.080760915236984d0*x490 - 0.0302853
     &43213869d0*x491) + x317*x461* x576*x580 + x318*x585 + 4.2666666666
     &6667d0*x331*x586 + 0.807691733333333d0*x574*x575 - x574*(-0.005090
     &39121236282d0*x169 + 0.00190889670463606d0*x172 + 0.00190889670463
     &606d0*x174)) + 0.733761610865473d0*x125*x206 - 0.80760915236984d0*
     &x131*x242 - x133*x231*x481 - x153*x444*x445 + 0.489174407243648d0*
     &x161*x71 + x166*x492*x89 - x168*x456 + x168*x460 - 0.0945455706654
     &444d0*x175 *x468*x469 + x179*(x180*x475 + 0.0757133580346725d0*x44
     &0 + 0.740308389672353d0*x472 - 0.40380457618492d0*x473 - 0.4038045
     &7618492d0*x474 + x476*x49 + 0.0757133580346725d0*x477) - x183*x470
     &*x471 + x183*x496*x498 - x185*x481 - x186*(-x199*x535 - x201*x540 
     &+ x204*x540 + x536*x538) - x217*(x508 - 0.366880805432736d0*x509 +
     & x510*x73 + x530*x73) - x224*x492*x82 + x229*(-x226*x456 + x226*x4
     &60 + 0.0121861463320552d0*x482 - 0.000566174082457926d0*x483 + 6.3
     &5065792816135d-6*x485 + 0.164512975482745d0*x486 + 0.0897343502633
     &155d0*x487 + 0.0168251906743717d0*x488 - 0.0897343502633155d0*x490
     & - 0.0336503813487433d0*x491) - 2.0d0*x232*x493 - x243*x505 + x262
     &*( -0.366880805432736d0*x161 - x211*x248*x287 + x213*x248*x540 + 0
     &.550321208149104d0*x214*x40 - 0.166666666666667d0*x244*x287*x54 - 
     &x257*x544*(x545 + x546 - x547 + x549 - x551 + x552 + x553 + x554) 
     &+ x285*x286*x544 + 0.166666666666667d0*x543) - x266*(x204* x561 - 
     &x211*x536*x537 + x557 - x558*x567 + x564*x86 - x565*x87) + x267*x5
     &05 + x274*(-1.21141372855476d0*x473 + 0.227140074104017d0* x478 + 
     &0.454280148208035d0*x506 + x507) - x292*( 0.0841259533718583d0*x12
     &8*x277 + 0.45860100679092d0*x161 - 0.0630944650288937d0*x2*x275*x5
     &71 + 0.416666666666667d0*x244*x259 - x255*x289*x571 - x257*x289*(-
     &x545 - x546 + x547 - x549 + x551 - x552 - x553 - x554) + x278*x510
     & + x278*x530 + x279*x284*( -0.189283395086681d0*gamma_aa*x527 - 0.
     &378566790173362d0*x172* x520 - 0.841259533718583d0*x511 - 1.850770
     &97418088d0*x514 + 1.0095114404623d0*x516 - 0.189283395086681d0*x52
     &3*x525 - 0.189283395086681d0*x548 + 1.0095114404623d0*x550) + x279
     &*x573*( -0.264996753121354d0*x172*x193 + 0.70665800832361d0*x187 +
     & 0.70665800832361d0*x190 - 0.264996753121354d0*x197)/x58**4 - 0.68
     &7901510186381d0*x509 - 0.416666666666667d0*x543) + x417*x48* x566*
     &(x201*x561 - x211*x535 - x38*x558*x559 + x38*x564*x78 + x557 - x56
     &5*x80) + 1.48061677934471d0*x429*x47*x81 + x431*x432 + x432* x441 
     &- x442*x444*x81 - x461*x500*(0.30285343213869d0*x174 - x442 + x479
     &) - 0.0472727853327222d0*x462*x464 + x462*x465*x466 - x465* x496*x
     &531 + x502*x84*(x200*x237 + x234) - x504*(x204*x237 + x236 ) + 2.0
     &d0*x541*x542 + 2.0d0*x555*x556) + 7.26848237132856d0*x123* x333)
      d2Ea(2) = -x120*(x121*(0.015286700226364d0*gamma_aa*x145*x408*x419
     & - 0.000984849694431712d0*gamma_aa*x164*x175*x597*x637 + x101*( 0.
     &000188724694152642d0*x221 + x226*x620 + x227 - 0.0448671751316578d
     &0*x605 + 0.0168251906743717d0*x623 - 2.38149672306051d-6*x669 - 0.
     &0168251906743717d0*x670 + 0.0448671751316578d0*x671 + 0.0168251906
     &743717d0*x672 - 0.0168251906743717d0*x674) - x118*(x109*(0.0687901
     &510186381d0*x1* x374*x587 + x323*(x195*x730 + x224*x356 + x357*x72
     &9 - x360*x729 - x615 + x658 - x713) + x324*x728 - x329*x717 + x376
     &*x589*x591 - x591*x727 + x594*x639*x694) + x132*x317*x724 - x133*x
     &307*x373 - 0.00190889670463606d0*x293*x640 + x297*(x294*x722 + x29
     &5 - 0.227140074104017d0*x616 + x693 + 0.60570686427738d0*x711 - 0.
     &227140074104017d0*x723) + 0.134601525394973d0*x301*x523*x577* x719
     & + 0.40380457618492d0*x301*x720*x721 - x308*(x196*x726 + x306 *x66
     &0 + x307 + x313*x654 + x313*x658 - 0.179468700526631d0*x656 + 0.06
     &73007626974867d0*x657 - 0.179468700526631d0*x661) - 0.067300762697
     &4867d0*x312*x649 - x315*(0.000169852224737378d0* x221 + x310*x620 
     &+ x311 - 0.040380457618492d0*x605 + 0.0151426716069345d0*x623 - 2.
     &14334705075446d-6*x669 - 0.0151426716069345d0*x670 + 0.04038045761
     &8492d0*x671 + 0.0151426716069345d0*x672 - 0.0151426716069345d0*x67
     &4) + x331* x386 + 2.13333333333333d0*x381*x586 + 0.122306106775024
     &d0*x575* x720 - x584*(x154*x722 + x293 - x566*x713 - x617 + x664 +
     & 0.40380457618492d0*x711) + x724*x725) - x127*x223 - x127*x605 - x
     &132*x155 + x132*x178*(x173*x660 + x175 + 0.0757133580346725d0* x65
     &5 - 0.20190228809246d0*x656 + 0.0757133580346725d0*x657 + 0.075713
     &3580346725d0*x659 - 0.20190228809246d0*x661 + 0.0757133580346725d0
     &*x663) - x132*x231*x665 + x132*x273*(x268* x660 + x270 + 0.2271400
     &74104017d0*x655 - 0.60570686427738d0*x661 + 0.227140074104017d0*x6
     &63 + x693) - x132*x617*x81 + x152*x420* x642 + x155*x619 - x163*x4
     &20*x643 - x163*x635 + x166*x422 + x167* x623 + x181*x619*x81 - x18
     &3*x497*x645*x646 - x183*x639*x651*x652 + 0.227140074104017d0*x183*
     &x640*x647 - x184*x665 + x193*x495*x645 - x216*x409*(0.006944444444
     &44444d0*x682 - 0.0625d0*x684 + x692* x73) - x232*x650 + x242*x368*
     &x668*x678 - x267*x44*x680*x77 - 0.227140074104017d0*x271*x531*x648
     & + x3*x370*x653 + x314*x653 + 0.0138888888888889d0*x32*x392*x636 -
     & x372*x653 - x410*(-x249*x707 + 0.0104166666666667d0*x682 - 0.0937
     &5d0*x684 - 0.00946416975433406d0*x703 + 0.0157736162572234d0*x705 
     &+ 0.825481812223657d0*x706) + x422*x620 - x425*(x203*(x537*x694/x6
     &8 + x717) - x423*x715 + x567*x680 + x676*x709 - x679*x710*x86) + x
     &44*x504*x676 - 0.00715836264238522d0*x464*x640 - x465*x599*x646* x
     &677 - 0.00357918132119261d0*x469*x648 - 0.151426716069345d0*x471 *
     &x649 - x493*x602 - x494*x666 + 0.151426716069345d0*x498*x648 - 0.3
     &0285343213869d0*x500*x640 - x502*x667*x668*x82 + x541*x702 + 0.227
     &140074104017d0*x542*x596 + x555*x708 + x556*x701 - x598*( x114*x39
     &6*(-x281*x662 - x405*x569 - x568 + x578 - x579 - 0.189283395086681
     &d0*x655 + 0.50475572023115d0*x656 - 0.189283395086681d0*x657 - 0.1
     &89283395086681d0*x659 + 0.50475572023115d0*x661) - x288*x707 - 0.0
     &0501594851177569d0*x395 *x408*x573*x718 + 0.00597136727592344d0*x3
     &96*x603 + 0.0130208333333333d0*x682 - 0.1171875d0*x684 - 0.0236604
     &243858352d0*x703 + 0.0394340406430586d0*x705 + 1.03185226527957d0*
     &x706) - x700*(-x201*x699 + x204*x699 - x536* x697 + x695*x696) + x
     &716*(-x2*x414*x710*x78 + x200*(x38*x714 + x534*x694/x65) + x201*x6
     &51*x709 + x414*x559*x651 - x60*x715*x80 )) + x124*(-x118*x595 - x2
     &31*x602 - x389*x408 + x390*x596 - x393* x44*x597 + x411*x600 + x41
     &2*x596 + x413*x601 - x420*x604 + x426 + x598*x599))
      d2Ea(3) = x122*(-0.011465025169773d0*gamma_aa*x356*x643 + 0.000111
     &849416287269d0* mu*x115*x53*x755 - x101*(0.0336503813487433d0*x739
     & - 0.0168251906743717d0*x749 - 0.0336503813487433d0*x750 - 0.03365
     &03813487433d0*x751 - 4.71811735381605d-5*x776 + 8.93061271147691d-
     &7*x777 + 0.0168251906743717d0*x778 + 0.0168251906743717d0*x779) + 
     &x118*(x109*(-0.0687901510186381d0* x105*x746 + x322*(-2.0d0*x356*x
     &730 + x736 - x795 + x798) + x375* x377 + x376**2*x589 - x376*x727 
     &+ x378*x728 - x379*x800 + 0.00143312814622163d0*x593*x774) + 0.035
     &6689671948494d0*x29*x577* x770 + x296*x364*(-0.454280148208035d0*x
     &357 + x780 - 0.227140074104017d0*x796 - 0.454280148208035d0*x797 +
     & 0.227140074104017d0*x799) + x299*x317*x364*x801 - x315*( 0.030285
     &343213869d0*x739 - 0.0151426716069345d0*x749 - 0.030285343213869d0
     &*x750 - 0.030285343213869d0*x751 - 4.24630561843445d-5*x776 + 8.03
     &755144032922d-7*x777 + 0.0151426716069345d0*x778 + 0.0151426716069
     &345d0*x779) - x364* x583*(-0.30285343213869d0*x357 + x767 - 0.1514
     &26716069345d0*x796 - 0.30285343213869d0*x797 + 0.151426716069345d0
     &*x799) - x365*( 0.000162998008060326d0*x3 + x313*x789 + 2.37463589
     &891975d-6*x343 - 1.70225887187107d-5*x346 + 8.10221651407408d-5*x3
     &49 + 0.000358795087982141d0*x352 - 0.000247481481481481d0*x368 + 0
     &.134601525394973d0*x406 + x726*x761 + 0.0673007626974867d0*x737 + 
     &0.134601525394973d0*x760) - 0.134601525394973d0*x373*x596 + 1.1306
     &5281331778d0*x381*x385 + x382*x725*x801 + 0.0277806181562528d0*x52
     &*x771 + 0.107006901584548d0*x721*x771 - 0.000289057959298633d0*x77
     &4) + x167*x749 - x178*x364*( 0.000183372759067866d0*x3 + 2.6714653
     &8628472d-6*x343 - 1.91504123085495d-5*x346 + 9.11499357833334d-5*x
     &349 + 0.000403644473979909d0*x352 - 0.000278416666666667d0*x368 + 
     &0.151426716069345d0*x406 + 0.0757133580346725d0*x737 + 0.075713358
     &0346725d0*x759 + 0.151426716069345d0*x760 + 0.0757133580346725d0*x
     &763) - x201*x501*x775 - x203*x503*x772 + x231*x364*x768 + x233*x41
     &4*x77*x775 + x266*x678*x772 + 0.454280148208035d0*x271*x400*x677*x
     &757 - x273*x364*( 0.454280148208035d0*x406 + 0.227140074104017d0*x
     &759 + 0.454280148208035d0*x760 + 0.227140074104017d0*x763 + x780) 
     &- x356*x635 + 0.011465025169773d0*x361*x642 - x364*x738*x81 + x388
     &* x768 - 0.0260152163368769d0*x409*x53*x766 - x410*( 0.82548181222
     &3657d0*x766 - 0.00354906365787527d0*x787 + 0.00179141018277703d0*x
     &788 - 0.0833333333333333d0*x790) + x425*( -x203*x60*x791*x794 + x2
     &12*x792*x794*x86 - x536*x793 + x675*x800 - x679*x800*x87 + x793*x8
     &7*x97) - x445*x739 + x445*x750 + x445* x751 + 0.000541983673684936
     &d0*x463*x756 + 0.30285343213869d0*x497 *x644*x757 + 0.068790151018
     &638d0*x499*x756 + 0.30285343213869d0* x596*x650 - 0.45428014820803
     &5d0*x596*x702 + x598*( 0.000949435763888889d0*rho_a**(-7.0d0)*x572
     &*x718*x781 - 0.0189283395086681d0*x397*(0.000458431897669666d0*x3 
     &+ 6.67866346571181d-6*x343 - 4.78760307713738d-5*x346 + 0.00022787
     &4839458333d0*x349 + 0.00100911118494977d0*x352 - 0.000696041666666
     &667d0*x368 + 0.378566790173362d0*x406 + 0.189283395086681d0*x737 +
     & 0.189283395086681d0*x759 + 0.378566790173362d0*x760 + 0.189283395
     &086681d0*x763) - 1.03185226527957d0*x766 + 0.00443632957234409d0*x
     &787 - 0.00447852545694258d0*x788 + 0.208333333333333d0*x790) - 2.0
     &d0* x600*x601*x666 + 0.011465025169773d0*x604 + 0.0001491325550496
     &92d0*x637*x755 - 0.034395075509319d0*x647*x756 + 0.011465025169773
     &d0*x652*x774 + x700*(-x30*x536*x782 + x675* x765 + x785) - 2.0d0*x
     &701*x708 + x716*(-x22*x784*x79*x97 - x29* x78*x786*x792 + x414*x66
     &7*x781*x791 + x415*x783 + x785))
      end subroutine
